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ABSTRACT 

Aims. (Ultra) Luminous Infrared Galaxies [(U)LIRGs] are important galaxy populations for studying galaxy evolution, 
and responsible for a significant fraction of the past star formation beyond z ~ 1. Local (U)LIRGs can be used to study 
criteria adequate to characterize similar high-z populations. Of particular interest is to identify reliable kinematic-based 
methods capable of distinguishing disks and mergers, as their relative fraction is a key observational input to constrain 
different evolutionary scenarios. 

Methods. Our goal is to analyze in detail the kinematics of the ionized gas of a small sample of LIRGs and study criteria 
able to characterize the status of these systems. 

Results. We have obtained VIMOS/VLT optical integral field spectroscopy (IFS) data for a sample of 4 LIRGs which 
have been selected at a similar distance (~ 70 Mpc) to avoid relative resolution effects. They have been classified in two 
groups {isolated disk and post-coalescence mergers) according to their morphology. The kinemetry method (developed 
by Krajnovic and coworkers) is used to characterize the kinematic properties of these galaxies and to discuss new 
criteria for distinguishing their status. 

Conclusions. We present and discuss new kinematic maps (i.e., velocity field and velocity dispersion) for these four 
galaxies. When the full 2-D kinematic information of the IFS data is analyzed with kinemetry, their morphological and 
kinematic classifications are consistent, with disks having lower kinematic asymmetries than post-coalescence mergers. 
We propose and discuss a new kinematic criterion to differentiate these two groups. In particular, we introduce a weight- 
ing which favors the outer parts of the kinematic maps when computing the total asymmetries. This is motivated by 
the fact that post-coalescence mergers show relatively small kinematic asymmetries in the inner parts as a consequence 
of the rapid relaxation into a rotating disk, with the outer parts being still out of equilibrium (i.e., larger asymmetries). 
This new criterion distinguishes better these two categories and has the advantage of being less sensitive to angular 
resolution effects. According to the previous kinemetry-based work aimed at distinguishing disks and mergers at high-z, 
the present post-coalescence systems would have been classified as disks. This indicates that the separation of disks 
from mergers is subjective to the definition of 'merger'. It also suggests that previous estimates of the merger/disk ratio 
could have been underestimated, but larger samples are necessary to establish a firmer conclusion. 

Key words, galaxies - kinematics - luminous infrared galaxies - integral field spectroscopy 



1. Introduction 

The kinematic characterization of high-z galaxy popula- 
tions is a key observational element to distinguish between 
different galaxy evolutionary scenarios, since it may help 
to find out the fraction of rotating disks and mergers at 
different cosmic epochs. This is a key observational input 
to constrain the relative role of major mergers vs. steady 
cool gas accretion in shaping galax i es, w hi ch still is un- 
der discussion (e.g. iGenzel et alj [200ll: iTacconi et al 
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With the aid of integral field spectroscopy (IFS) we can 
resolve the kinematic status and internal processes at work 
within galaxies and better understand the role of the dom- 
inant mechanisms involved at early epochs. For instance, 
iForster S chreibe r et aD (|2009[ ) classified the SINS sample 
(i.e., 63 galaxies with 1.3 < z < 2.6) concluding that these 
objects can be classified as: 'rotation-dominated', interact- 
ing or merging systems and 'dispersion-dominated' (i.e.. 



the majority of the star forming galaxies at these redshifts 
is fed by gas via continuous cold flows along streams (in- 
cluding minor mergers). 

At intermediate redshift (z ^ 0.6), spatially resolved 
kinematics has revealed a large fraction of chaotic veloc- 
ity fields, supporting that most loca l spirals were rebuilt 
after a maj or merger since z = 1 (iHammer et al. 20051 : 
Puechll20ldf). At s imilar redshifts, iFlores et al.l ()2006h and 
Yang et aLl (|2008D investigated the velocity fields of galax- 
ies with IFS, and derived that only 35% of their sources can 
be classified as rotating disks while the rest show perturbed 
or complex kinematics. 

Some local galaxy populations are particularly rele- 
vant for the study of galaxy evolution. This is the case 
of nearby Luminous and Ultraluminous Infrared Galaxies 
(LIRGs, Lib. = [8-1000 ^lm] =10" - lO^^ Lq, and ULIRGs, 
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Lju > 10^^ Lq, respectively). Despite these objects are 
relatively rare in the local universe, they are much more 
numerous at high-z and are responsible of a significant frac- 
tion of the whol e pas t star formation beyon d z ~ 1 (e.g ., 
iLe Floc'h etlD [2005t iPerez-Gonzalez et all [20051 mm . 
These objects exhibit a large range of interacting/merging 
properties, from mostly isol ated discs for low l uminosity 
LIRGs (e.g. log LiR < 11.3. lArribas et al.ll2004) to a ma- 
jority of merger rem nants for ULIRGs (e.g. log Lju > 12, 
IVeilleux et al.ll200a) . Local (U)LIRGs were initially consid- 
ered as the local counterpart of high-z (U)LIRGs like those 
discovered by Spitzer and the more lu minous sub millimiter 
galaxies (SMGs; e.g. Blain et al. 200llTacconi et al. 20061. 
The study of local (U)LIRGs is also important to better 
understand the star formation history in the universe and 
it can provide an opportunity to correlate the properties 
of local and high-z populations. Moreover, large scale tur- 
bulences, tidal tails and outflows can be studied in detail 
allowing a good characterization of these objects. However, 
there are some discussion about the similarities between lo- 
cal and distant (U)LIRGs (e.g.. lRuiopakarn et al.ll2011l and 
references therein). 

Recently several authors have analyzed the velocity 
fields and velocity dispersion maps of different galaxy 
samples with the aim of discerning merging and non- 
merging systems based on their k inem atic pr operties us - 
ing the kinemetry methodol ogv (IKra inovic e t al.l 1200611 . 
iJesseit et al.l ()2007D and iKronberger et al.i (i2007l ) have in- 
vestigated the distortions in the velocity fields of simulated 
interacting disc/merger galaxies at differen t redshifts in 
the range z — — 1, while Shapiro et alj (|2008) (here- 
after, SOS) provide a straightforward mean of classifying 
z ^ 2 galaxies, observed with SINFONI/VLT, as disks or 
mergers based on the asymmetries of their velocity field 
and velocity dispersion maps. They compare and distin- 
guish two kinematic classes characterized by recent major 
merger events (i.e., mergers) and those without signs of in- 
teracting or recent merger activity (i.e., disks). 

Similar work s based on IFS data (i.e ., OSIRIS/Keck, 
iBasu-Zvch eTall l2009t iGonpalves etall l201dl) are fo- 
cused on nearby objects (i.e., supercompact UV-Luminous 
Galaxies (ScUVLGs) and Lyman Break Analogs (LBAs)). 
The morphology and kinematics of LBAs projected at high 
redshift show similarities w ith compact and disp ersion- 
dominated z ^ 2 galaxies (jBasu-Zvch et al.l l2009h . This 
suggests that galaxy interaction and/or mergers could also 
be driving the dynamics of dispersion-dominated z '^ 2 
galaxies. 

In this paper we present new spatially resolved kine- 
matics of four local LIRGs obtained with the VLT/VIMOS 
IFU as part of a larger project to characterize the properties 
of (U)LIRGs on th e basis of optical and infrared Integral 
Field Spectroscopv (ICohna et a l. 2005: Arri bas et al.l[200l 
iMonreal-Ibero et al.l 120101: iRodriguez-Zaurin et al.l 120111 
and references therein). The fact that two of our LIRGs 
are post-coalescence mergers will offer the opportunity to 
explore the potential of the kinemetry method when an- 
alyzing the velocity fields and velocity dispersion maps in 
moderately disturbed and partially relaxed systems. The 
paper is structured as follows. In Section 2, we present the 
sample selection as well as some details about the obser- 
vations, data reductions, line fitting and velocity map con- 
struction. Section 3 is devoted to the description of the main 
kinematic properties of the individual objects. In Section 



4, the kinemetry analysis and its potential to distinguish 
discs and mergers are presented and discussed with the 
aid of different asymmetry planes. Finally, the main results 
and conclusions are summarized in Section 5. 

Throughout the paper we will consider Hq = 70 
km/s/Mpc, riM = 0.3 and flA = 0.7. 



2. The VIMOS sample, observations, data 
reduction, line fitting and map construction 

2.1. The sample and morphological classification 

For the present analysis we have selected 4 LIRG galax- 
ies at a similar distance (^ 70 Mpc): two of them (i.e., 
IRAS F11255-4120, IRAS F10567-4310) are morpholog- 
ically classified as class objects and two (i.e., IRAS 
F04315 -0840, IRAS F 2 1453-3 511) as class 2 objects accord- 
i ng to lArribas et al.l (|2008[ ) and iRodriguez Zaurin et al.l 
()2010[ ) (hereafter. Paper I and Paper III respectively). In 
this classification scheme class is defined as an object 
with relatively symmetric morphology, appearing to be iso- 
lated with no evidence for strong past or ongoing interac- 
tion (hereafter, disks), while class 2 is an object with a 
morphology suggesting a post-coalescence merging phase, 
with a single nucleus or two distinct nuclei at a projected 
distance D < 1.5 kpc (hereafter, post-coalescence mergers). 
We have chosen these four systems at the same distance as 
they will allow us to discuss their relative kinematic proper- 
ties avoiding linear resolution dependency effects. In tablejT] 
we present the main properties of the sample: note that the 
systems classified as post-coalescence mergers have higher 
IjiR than the disks. 



2.2. Observations 

The observations were carried out in service mode usin g 
the Integral Field Unit of VIMOS (|Le Fevre et al.ll2003h . 
at the Very Large Telescope (VLT), covering the spectral 
range (5250 — 7400) A with the high resolution mode 'HR- 
orange' (grating GG435) and a mean spectral resolution of 
3470. The field of view (FoV) in this configuration is 27" 
X 27", with a spaxel scale of 0.67" per fiber (i.e., 1600 
spectra are obtained simultaneously from 40 x 40 fibers 
array). A square 4 pointing dithering pattern is used, with 
a relative offset of 2.7" (i.e., 4 spaxels). The exposure time 
per pointing ranges from 720 to 850 seconds, such that, 
the total integration time per galaxy is between 2880 and 
3400 seconds. For further details about the observations see 
Paper I. 

2.3. Data reduction 

The VIMOS data are reduced with a combination of the 
pipeline Esorex (version 3.5.1 and 3.6.5) included in the 
pipeline provided by ESO, and different customized IDL 
and IRAF scripts. The basic data reduction is performed 
using the Esorex pipeline (bias subtraction, fiat field cor- 
rection, spectra tracing and extraction, correction of fiber 
and pixel transmission and relative flux calibration). Then 
four quadrants per pointing are reduced individually and 
combined in a single data cube associated for each point- 
ing. After that, the four independent dithered pointings are 
combined together to end up with the final 'super-cube'. 
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Table 1. General properties of the LIRGs sub-sample. 



IDl 


ID2 


a 


.5 


^ 


D 


scale 


log L/H 


Class 


References 


IRAS 


Other 


(J2000) 


(J2000) 




(Mpc) 


(pc/arcscc) 


(Lo) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


F11255-4120 


ESO 319-G022 


11:27:54.1 


-41:36:52 


0.016351 


70.9 


333 


11.04 





1 


F10567-4310 


ESO 264-G057 


10:59:01.8 


-43:26:26 


0.017199 


74.6 


350 


11.07 





1 


F04315-0840 


NGC 1614 


04:33:59.8 


-08:34:44 


0.015983 


69.1 


325 


11.69 


2 


1, 2 


F21453-3511 


NOG 7130 


21:48:19.5 


-34:57:05 


0.016151 


70.0 


329 


11.41 


2 


1, 2 



Notes. Col (1): object designation in the IRAS Faint source catalogue (FSC). Col (2): other name. Col (3) and (4): 
right ascension (hours, minutes and seconds) and declination (degrees, arcminutes and arcseconds) from the IRAS FSC. 
Col (5): redshift of the IRAS sources from the NASA Extragalactic Database (NED). Col (6): luminosity distances 
assuming a ADCM cosmology with Ho — 70 km/s/Mpc, Qm = 0.3 and Ha ~ 0.7, using the Edward L. Wright Cosmology 
calculator. Col (7): scale. Col (8): infrared luminosity (L/_r= L( 8-1000) /xrn), i n unit s of solar bolometric luminosity, 
calculated using the flu xes in the four IR AS bands as given in ISanders et al.l (|2003l ) when available. Otherwise the 
standard prescription in iSanders fc Mirabel (1996 ) with the values in t he IRAS Point and Fain t sourc e catalogues was 
used. Col (9): Morphology class. For further information see table 1 in [Rodriguez Zaurin et ahl l|2010ll . 

References. (1) iVeiUeux et all (|1995D : (2) iCorbett et all (|2003l ). 



containing 44 x 44 spaxels for each object (i.e., 1936 spec- 
tra). Further details about the data reduction for the whole 
sample can be found in Paper I and in Paper III. 

The wavelength calibration, the instrumental profile 
and fiber-to-fiber transmission correction are checked for 
the four galaxies using the [OIJA6300.3 A sky line as 
in Paper I. The results is focused on the Ha and 
[NIIJAA6548.1, 6583.4 A emission lines, so that the sky line 
here considered is suitable for its proximity to these lines. 
In table [H the average values for the central wavelength 
< Ac > and (instrumental) width < ajNS >: with their 
standard deviations, are shown for the four galaxies. 



Table 2. Results of the calibration check using the 
[OI]A6300.3 A sky line. 



IRAS Galaxy 


< A, > ±AA, [A] 


< UINS > ±std [A] 


F11255-4120 


6300.42 ± 0.15 


0.78 ± 0.09 


F10567-4310 


6300.17 ± 0.13 


0.77 ± 0.09 


F04315-0840 


6300.24 ± 0.10 


0.78 ± 0.08 


F21453-3511 


6300.27 ± 0.11 


0.77 ± 0.06 



Notes. Typical values of the central wavelength and width dis- 
tribution of the [OI]A6300.3 A sky line registered in the data 
cubes along with their standard deviations (errors). 



of the fit. For each emission line we end up with the fol- 
lowing information: central wavelength (Ac), intrinsic width 

(aune) (i.e., aiine = V'^lbs - ^Ins) ^nd flux intensity along 
with their respective errors. To obtain the radial velocity 
uncertainty estimates, wavelength calibration and fitting 
errors were combined in quadrature, giving a global wave- 
length error of AAtot 



^aaT 



''MP FIT + ^^c- For high S/N 
spectra the wavelength calibration error can be considered 
the main source of uncertainty since the fitting errors are 
typically substantially smaller (and vice versa for low S/N 
spectra). For the fitting errors, we derive values generally 
smaller than 0.2 A while the wavelength calibration errors 
are of the order of ^ 0.12 A as shown in Tab. [31 An average 
value of the ajNS of (85 ± 9) km/s is derived at the corre- 
sponding wavelength, along with a mean FWHM of (1.82 
± 0.19) A . 

We start fitting all the lines to single Gaussian profiles 
since the main component usually extends over the entire 
galaxy. Two-Gaussian profiles (i.e., main and broad compo- 
nents) are needed for the inner regions of the four objects. 
Using a c ollection of procedur es written in IDL code (i.e, 
imaplot. lMaiz-ADellanial2004l ) ■ we generate flux intensity, 
velocity field and velocity dispersion maps, respectively, for 
the main and broad components (Figures [T]- 2]) along with 
a VIMOS continuum image. When HST images are avail- 
able (i.e, NICMOS F160W (H-band) and ACS F814W (I 
band)) they are also shown in the panels. 



2.4. Line fitting and map construction 

As mentioned in the previous paragraph, we are interested 
in the Ha - [Nil] emission line complex analysis. The ob- 
served emission profiles of the individual spectra are fit- 
ted to a Gaussian function using an IDL routine (i.e., 
MPFITEXPR, implemented by C. B. Markwardt). This al- 
gorithm finds the best set of model parameters which match 
the data and it is able to fix the line intensity ratios and 
the wavelength differences according to the atomic param- 
eters when adjusting multiple emission lines (i.e., the Ha - 
[Nil] complex). The MPFIT routine provides errors for the 
output parameters and gives an estimate of the goodness 



3. Kinematic properties of individual objects 

In this section we describe the global kinematic behavior as 
inferred from our IFS maps (i.e, fiux intensity, velocity field 
and velocity dispersion). We will also infer some kinematic 
parameters (e.g., ac, Vc/o'c, ^shear/^ and dynamical mass 
M.dyn) useful to characterize these systems. 

The ratio of the maximum circular rotation velocity Vc 
(i.e., the half of the observed peak-to-peak velocity) and 
the central velocity dispersion CTc measures the nature of 
the gravitational support of a system in equilibrium. A 
"^cl <^c > 1 is the signature of a rotation-dominated sys- 
tem while a lower value (i.e., \cl<^c < 1) means that the 
object is dispersion-dominated, as in the case of elliptical 
galaxies, where the galaxy is sustained against gravitational 



Bellocchi et al.: Kinematic asymmetries of disks and post-coalescence mergers 



collapse by the pressure originated by random motions of 
the stars. When analyzing the ^d <yc parameter it is im- 
portant to realize that the presence of flows, superwinds or 
other merger-induced processes are factors that may lead to 
large Vc and/or Oc values. Therefore we compute the circu- 
lar velocity Vc and the central velocity dispersion CTc using 
the component of one-Gaussian model fit (i.e., Ic) and the 
main (i.e., the one defining the systemic behaviour) com- 
ponent of the two-Gaussian model fit (i.e., 2c). 

Since class objects show the features of an ideal ro- 
tating disfcQ we correct their kinematic values for their re- 
spective inclination; class 2 objects show the velocity fields 
distorted such that we derive their kinematic values with 
and without correcting for the inclination. 

The ratio between the velocity shear y shear and the 
global velocity dispersi on in the whole galaxy S has been 
derived too, as done in lGoncalves et all ()201Cl[ ). The ^ shear 
(not corrected for the inclination of the object) has been 
computed as the median of the 5-percentile at each end of 
the velocity distribution, for the Vmax and v„ii„ and then 

QenneQ as V shear — 2 V^"^*i2; ^min)' 

We have also used the kinematic information to derive 
the dynamical mass M^yn- This mass takes into account the 
whole gravitational field and so dark matter, gas and stel- 
lar components are included. Spiral galaxies can be modeled 
using a fiat component for the disk and a spheroidal com- 
ponent for the central concentration (i.e., a de Vaucouleurs 
profile and a massive halo). Assuming that the central re- 
gions are virialized, the observed velocity dispersion of the 
ionized gas Uc can provide a good estiniate fo r the disk's 
dynamical mass {Mdyn, e.g.. lColina et al.l[2005[ ). We use the 
following relation: 



M, 



dyn 



m W Rhm < Mq 



(1) 



3.1. IRAS F11255-4120 (ESO 319-G022) 

This is a class object, appearing to be a single isolated 
galaxy with a relatively symmetric morphology according 
to its DSS image. It is a barred spiral with a ring extended 
up to ~ 4 kpc from the nuclear region, clearly detected in 
our Ha image. The orientation of the bar is different in 
the VIMOS continuum image (P.A. ~ 110°) and in the Ha 
image (P.A. ~ 150°). 

For the main component, as shown in the top panel 
of Fig. m the velocity dispersion (a) map has a (almost) 
centrally peaked structure with values of 120 km/s (high 
values also in the bar structure) and a lower mean value (~ 
20 km/s) is found in the ring. The location of this peak does 
not agree with those of the Ha flux intensity and continuum 
maps with an offset of about 0.6 kpc (i.e., ~ 2 arcsec). The 
kinematic center of the velocity field well agrees with the 
Ha fiux peak. The projection of the rotation axis (minor 
kinematic axis) seems to be shifted by about 30° with re- 
spect to the orientation of the bar in the Ha fiux intensity 
map and in agreement with the bar of the continuum. The 
observed velocity amplitude for the main component at a 
distance of 4 kpc from the H^ peak is of the order of (300 
± 8) km/s. 

A small region of about ~ 0.67 kpc x 0.89 kpc in the 
nuclear region shows traces of a second component, indi- 
cating the presence of a non-rotational mechanisms (e.g., 
outflow / wind). This component is blueshifted by 80 km/s 
with respect to the main one and shows a broad velocity 
dispersion a in the range (160 - 220) km/s and a velocity 
amplitude of the order of (60 ± 7) km/s. 

The yc/<7c and Vshear/^ have values of 2.3 and 2.7 re- 
spectively, supporting the idea that this galaxy is rotation- 
dominated: the kinematic properties of this galaxy seem to 
be consistent with its morphology (i.e., a rotating disk). 



where Rhm is the half-mass radius in kpc and CTc is the 
central velocity dispersion in km/s. The Rhm parameter is 
not an observable itself and therefore it cannot be measured 
directly, so it can be inferred from measuring the radius 
that encloses half luminosity (i.e., Reff)- The Re// were 
obtained from existing 2MASS near-infrared imaging. The 
extinction in the near-IR is a factor of 5-10 smaller than 
that in the optical and a proper way to minimize these ef- 
fects is to measure the half-light radius in the H ot K bands 
as done here. The parameter m depends on the mass dis- 
tribution: its value ranges from 1.4, for a King stellar mass 
distribution with a tidal-to-core rat io of 50, which is a g ood 
representation of elliptical galaxies (JTacconi et al.ll2002n . to 
1.75 for a poly tropic sphere with a density index covering a 
range of values (S pitzcr 1987), to 2.1 f or a de Vaucouleurs 
mass distribution (jCombes et al.lll995l ). As in lColina et al.l 
()2005() we will assume m = 1.75. All these results are shown 
in the table [31 their inclinations are drawn from the NED 
site and are consistent with our own estimates from the Ha 
distribution. 



^ We refer to ideal rotating disk as a thin disk with gas clouds 
kinematically characterized by: a velocity field that peaks at the 
galaxy major axis and goes to zero along the minor axis of each 
orbit, while the velocity dispersion will be constant long each 
orbit and will decrease between orbits of increasing major axes 
(see description in Sec. 4.1) 



3.2. IRAS F10567-4310 (ESO 264-G057) 

From the DSS image IRAS F10567-4310 is morphologically 
classified as class object (see Paper HI) . The Ha flux map 
inferred from our IFU data (top panel of Fig.[2|) shows a ring 
with a patchy distribution and its peak is in good positional 
agreement with the center of the of galaxy as inferred from 
the continuum. 

The main component shows a centrally peaked velocity 
dispersion map with typical values of 60 - 70 km/s. Its gen- 
eral morphology resembles that of the Ha flux map (i.e., 
regions with high Ha surface brightness tend to have large 
velocity dispersion, with a mean value of 40 km/s). This 
component shows a disk-like regular velocity field with a 
clear rotation pattern, with the kinematic center in a good 
positional agreement with the Ha fiux peak. The observed 
velocity amplitude is of (300 ± 6) km/s computed at a dis- 
tance of 4.7 kpc from the kinematic center. 

A secondary kinematic distinct component is identified 
in the nuclear region over an extension of ~ 1.2 kpc x 0.9 
kpc; it shows a velocity amplitude of (128 ± 9) km/s and it 
is blue-shifted by ~ 70 km/s with respect to the main one; 
its velocity dispersion is in the range (90 - 180) km/s. 

The Vc/o'c parameter clearly classifies this object as 
rotation-dominated with Vc/o'c '^ 4.6 and Vshear/^ ^ 3. 
Kinematics of this class object is also consistent with its 
morphology. 
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Fig. 1. Top panel: Ha maps for the class object IRAS 
F11255-4120. VIMOS continuum image; HST continuum images 
are not available for this galaxy. Middle row: the flux intensity, 
velocity dispersion a (km/s) and velocity field v maps (km/s) 
for the main component. On the bottom: respective maps for the 
broad component. Note the different FoV between the main and 
the broad components. The latter has been zoomed in since it 
covers a small area. The flux intensity maps are represented in 
logarithmic scale and in arbitrary flux units. All the images are 
centered using the Hq flux intensity peak and the iso-countors 
of the Ha flux are overplotted. Bottom panel: Ha- [Nil] ob- 
served spectra of IRAS F11255-4120 for selected inner regions 
(indicated by the coordinates in the top label using the same 
reference system as in the top panel) where the main and broad 
component coexist. The red curve shows the total Ha- [Nil] com- 
ponents as obtained from multi-components Gaussian flts. The 
green and blue curves represent respectively the main and broad 
components. 



Fig. 2. Top panel: Ha maps for the class object IRAS 
F10567-4310. VIMOS continuum image; HST continuum images 
are not available for this galaxy. Middle row: the flux intensity, 
velocity dispersion a (km/s) and velocity fleld v maps (km/s) 
for the main component. On the bottom: respective maps for the 
broad component. Note the different FoV between the main and 
the broad components. The latter has been zoomed in since it 
covers a small area. The flux intensity maps are represented in 
logarithmic scale and in arbitrary flux units. All the images are 
centered using the Hq flux intensity peak and the iso-countors 
of the Hq flux are overplotted. Bottom panel: Ha-[NII] ob- 
served spectra of IRAS F10567-4310 for selected inner regions 
(indicated by the coordinates in the top label using the same 
reference system as in the top panel) where the main and broad 
component coexist. The red curve shows the total Ha- [Nil] com- 
ponents as obtained from multi-components Gaussian fits. The 
green and blue curves represent respectively the main and broad 
components. 



3.3. IRAS F04315-0840 (NGC 1614) 

This is a class 2 object (Fig. [21) with a relatively asymmetric 
morphology: the DSS image shows a tidal tail extending for 
13 kpc from the nuclear region. This is a well studied, post- 
coalescence late merger, with a bright, spiral structure at 



scale of few kpc (1-3 kpc). The spiral structure in Ha 
flux map shows a different orientation with respect to the 
continuum. 

The velocity dispersion map of the main component has 
an irregular structure. Its peak (i.e., 220 km/s) is found at 
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Fig. 3. Top panel: Ha maps for the class 2 object IRAS 
F04315-0840. VIMOS continuum image and HST continuum im- 
ages (i.e., H and I bands). Middle row: flux intensity, velocity 
dispersion a (km/s) and velocity field v maps (km/s) for the 
main component. On the bottom: respective maps for the broad 
component. Note the different FoV between the main and the 
broad components. The latter has been zoomed in since it cov- 
ers a small area. The flux intensity maps are represented in 
logarithmic scale and in arbitrary flux units. All the images are 
centered using the Ha peak and the iso-countors of the Ha flux 
are overplotted. Bottom panel: Hq-[NII] observed spectra of 
IRAS F04315-0840 for selected inner regions (indicated in the 
top label using the same reference system as in the top panel) 
where the main and broad component coexist. The red curve 
shows the total Ha- [Nil] components as obtained from multi- 
components Gaussian flts. The green and blue curves represent 
respectively the main and broad components. 



necting the N-E to S-W part) is not aligned with that of 
the inner part (aligned as N-S direetion). 

A broad component is found in the inner region and 
covers a relatively large area of ^ 2.4 kpc x 2.7 kpc. This 
component is blue-shifted by 50 km/s with respect to the 
main component. Its velocity dispersion is in the range (70 
- 400) km/s with a velocity amplitude of (760 ± 13) km/s. 
The fact that the projection of the kinematic axis of this 
broad component has a shift of almost 90° with respect to 
the main component and the blue-shifted region of the ve- 
locity field shows the largest velocity dispersion (i.e., '^ 400 
km/s) supports the hypothesis of a dusty outflow, where 
the receding components (which are behind the disk) are 
obscured, making the whole profile relatively narrow with 
respect to the approaching component. 

We derive a vd (Jc > 1 even when the inclination cor- 
rection is not included (i.e., Vc/cfc ~ 1-4 - 2.3): it shows the 
dominance of an intrinsic rotation over random motions. 
Therefore, this parameter would classify this object as ro- 
tation dominated. 



3.4. IRAS F21453-3511 (NGC 7130) 

IRAS F21453-3511 is a peculiar class 2 object, with traces 
of spiral and asymmetric morphology from the HST and 
DSS images. The ionized gas, as traced by the Ha emis- 
sion, is concentrated in the nuclear region and mostly in 
the northern spiral arm. 

The asymmetric velocity dispersion map of the main 
component shows higher values in the northern arm and its 
central part (i.e., two local maxima can be revealed) with 
values of 80 km/s. The velocity field of the same component 
shows an asymmetric behavior where three main regions 
can be identified: the SW part having a mean value of 4900 
km/s, the NW one with 4850 km/s and the E part with 
4780 km/s. The NW part could have different inclination 
than the rest of the observed regions, so explaining this 
peculiar velocity structure. The kinematic center can be 
well identified with the Hq, flux peak although the rotational 
axis it is not well defined. The photometric major axis in the 
NICMOS image seems to be orientated in the N-S direction 
and the ACS image seems to reveal a ring of knots with a 
major axis orientated along P. A. ^ 135°. 

A small region of about 1.5 kpc x 1.3 kpc in the nuclear 
part shows a second component. The blue-shifted region 
of its velocity field shows the largest velocity dispersion 
(i.e., a ^ 400 km/s, see Fig. H]). This feature suggests that 
an outflow is present in the inner part of the galaxy. The 
amplitude of the velocity dispersion map is about (285 ±20) 
km/s while the velocity field amplitude is of (270±16) km/s. 
This component is blue-shifted by 150 km/s with respect 
to the main component. 

The derived Vc/ac and y shear /^ parameters are respec- 
tively 3 and 2: in this object the random motions do not 
seem to dominate although in its velocity field and veloc- 
ity dispersion maps there are some anomalies. With such a 
parameter we would classify this source as rotation domi- 
nated. 



2.4 kpc from the nucleus (i.e.. Ha flux peak) in the northern 
arm. The velocity field of the main component is somewhat 
distorted and chaotic with an amplitude of (325 ± 5) km/s. 
The projection of the rotation axis in the outer part (con- 



3.5. Summary of the global kinematics properties 

From the kinematic maps and the results derived so far 
(see table [3]) we can draw some general conclusions. Class 
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Fig. 4. Top panel: Ha maps for the class 2 object IRAS 
F21453-3511. VIMOS continuum image and HST continuum im- 
ages (i.e., H and I bands). Middle row: flux intensity, velocity 
dispersion a (km/s) and velocity field v maps (km/s) for the 
main component. On the bottom: respective maps for the broad 
component. Note the different FoV between the main and the 
broad components. The latter has been zoomed in since it cov- 
ers a small area. The flux intensity maps are represented in 
logarithmic scale and in arbitrary flux units. All the images are 
centered using the Ha peak and the iso-countors of the Ha flux 
are overplotted. Bottom panel: Hq-[NII] observed spectra of 
IRAS F21453-3511 for a selected inner region (indicated by the 
coordinates in the top label using the same reference systems 
as in the top figure) where main and broad component coexist. 
The red curve show the total Ha- [Nil] components as obtained 
from multi-components Gaussian fits. The green and blue curves 
represent respectively the main and broad components. 



objects are characterized by point-antisymmetric velocity 
fields, showing the spider diagram typical of a ideal rotat- 
ing disk; their velocity dispersion maps are centrally peaked 
with a radial decay at larger radii. On the other hand, 
class 2 objects show velocity fields which do not follow 



this pattern. Despite they arc to first approximation point- 
antisymmetric, they also show clear disturbances and asym- 
metries. These could be explained by a bent disk and/or 
other tidal motions, consequence of the past merger. Their 
a maps are more asymmetric with off-nuclear regions of pe- 
cuhar high velocity dispersion values (e.g., IRAS F04315- 
0840), and showing asymmetrical structures in the nuclear 
regions (e.g., IRAS F21453-3511). 

All of our sources show signs of a broad and blue-shifted 
components (i.e., with velocity shift of Av = (50 - 150) 
km/s with respect to the main component) in the inner 
part (nucleus) of the galaxy. They show mean velocity dis- 
persions ranging between (140 - 260) km/s, which are 2 - 
4 times higher than those of the main components. These 
components cover a small area in disk objects (i.e., ~ 1 
kpc^) while larger areas (i.e., 2-6 kpc^) are involved for 
the post-coalescence mergers, showing more complex and 
distorted spectra. Although this component is more promi- 
nent for the class 2 objects, the fact that it is blueshifted 
for the four sources and have high associated velocities 
suggests that an outflow can exist in their inner regions. 
Furthermore, the velocity field patterns of this component 
(i.e., kinematic axes perpendicular to those of the main 
component) support the dusty outflow hypothesis. 

From our kinematic results we see that the sole use of 
ID parameters Vcjoc and / or Yshear/'^ does not seem to 
give us a unique classification for our sources. Indeed, ac- 
cording to them, our objects seem to be rotation dominated, 
even when the inclination correction is not included. Their 
^c/o'c values are between thos e typical of local spiral galax- 
ies (i.e., Vc/ctc ~ 5 - 15. Epinat et al.ll2010) and those ob- 
tained f or Lyman Break Analo gs at z ~ 0.2 (with Vc/fXc ^ 
0.4 - l.S. lGonqalves et al.ll2010D . The derived global velocity 
dispersions E generally are bet ween (40 - 60) km/s, q uite 
comparable to those obtained in lGonpalves et all ()201 0^ for 
LBAs, with a typical median value of '^ 67 km/s, but much 
higher than those observed in other lower luminosity lo- 
cal star-forming ga laxies (i.e., typical u ~ 5 - 15 km/s, 
e.g.. lDib et al.ll2006l ). For a comparison we also derived the 
^ shear (not Corrected for the inclination of the galaxy) as 
explained before. Our values ranges between 50 (i.e., IRAS 
F04315-0840) and 130 km/s. This parameter reveals higher 
degree of rotation if compared with the values obtained for 
LBAs (i.e., ^ shear < 70 km/s) where in many cases a sig- 
nificant velocity gradient has not been observed, and no 
actual rotation can be identified. Then, we derive Vshear/"^ 
showin g values typical of rotation dominated objects (i.e., 
~ 2 — 3.lEpi nat et ani2010t ) with the exception of the galaxy 
IRAS F04315-0840 (see Tab.[3l). 

The derived dynamical masses classify our LIRGs as 
moderate mass systems, characterized by a mean value of 
(1.8 ± 0.4) • 10^° Mq. Our mass estimates agree with 
th ose obtained for LIRGs in iHinz fc RieQ (|2006l ) and 
iVaisanen et all ((2008) , which are in the range (10^° - 10") 

Mq. 

The post-coalescence mergers do not show extreme kine- 
matic asymmetries in their maps, making them good exam- 
ples to test the potential of kinemetry for characterizing 
moderate asymmetries. 

4. Kinemetry analysis 

The goal of this section is to investi gate the potentia l of the 
kinemetry method, developed by iKrainovic et al.l (|2006l) 
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(hereafter K06), when analyzing the kinematic maps of 
(U)LIRGs. In particular we want to investigate how power- 
ful this methodology is for studying the kinematic asymme- 
tries. First we will simply apply the kinemetry method to 
the sample, drawing preliminary conclusions on the kine- 
matics of these objects, th en we will app l y the same cri- 
teria as those proposed by IShaoiro et al.l 12008") and also 
explore the potential of a new criterion to study kinematic 
asymmetries. Therefore, we expect to find out lower asym- 
metries for our two post-coalescence merger systems with 
respect to those considered by SOS, which are more kine- 
matically disturbed due to recent merger activity. In our 
post-coalescence systems, the inner parts are expected to 
be relaxed into an almost virialized disk with a large rota- 
tional component while the outer parts should still retain 
asymmetries associated to the merger events. Finally we 
analyze the resolution / rcdshift effects on these results. 



of the ellipse {fJ.even{'>P,r) w ^o('')) and decreases when the 
semi-major axis length increases. In this context, the Aq 
term represents the velocity dispersion profile. Higher or- 
der terms (An, Bn) wi ll identify deviation from symmetry 
(jKrainovic et al.ll2006[ ). 

There are two effects which limit the reliability with 
which coefficients in the expansion can be determined: (i) 
the absolute number of points sampled along the ellipse, 
and (ii) the regularity with which these points sample the 
ellipse as a function of angle, 9. Ellipse parameters (i.e., 
position angle F, centre and flattening q — b/a) can be de- 
termined by minimizing a small number of harmonic terms; 
in the code some of them can be constrained (e.g., fix the 
center allows the position angle F and flattening q to freely 
change or remain fixed). The choice of a correct center is 
important in the analysis in order to avoid artificial over- 
estimations of the asymmetries. 



4.1. The method 



The k inemetry method, developed by iKrainovic et al.l 
(|2006f) . comprises a decomposition of the moment maps 
into Fourier components using ellipses. For clarity, we will 
briefiy describe again the main steps presented in K06 for a 
better understanding of this analysis. The Fourier analysis 
is the most straightforward approach to characterize any 
periodic phenomenon: the periodicity of a kinematic mo- 
ment can easily be seen by expressing the moment in polar 
coordinates: K (x, y) — )► K (r, ip). The map K(r, i/j) can be 
expanded as follows to a finite number (N+1) of harmonic 
terms (frequencies): 



N 



K{r,ij) ^Aoir) + ^Anir) sin{n-ij)+Bn{r) cos(n-V'), (2) 



ra=l 



where ip is the azimuthal angle in the plane of the galaxy 
(measured from the major axis) and r is the radius of a 
generic the ellipse. The harmonic series can be presented in 
a more compact way 



N 



K{r,^) = Aoir) + ^ fc„(r) • cos[n{^ - </)„(r))]. 



(3) 



where the amplitude and the phase coefficients (fc„, 
(pn) are easily calculated from the An, Bn coefficients: 

kn = y/A'^ + B'^ and 0„ = arctan ( -g^ ) . Thus, for an 

ideal rotating disk one would expect the velocity profile to 
be dominated by the Bi term while the velocity dispersion 
profile dominated by the Aq term. In the case of an odd 
moment (fJ-odd), the sampling ellipse parameters are deter- 
mined by requiring that the profile along the ellipse is well 
described by iioddiipjr) « Bi(r) ■ cos{ip) since the velocity 
field peaks at the galaxy major axis {tp — 0) and goes to 
zero along the minor axis (ip = 7r/2). So, the power in the 
Biy term therefore represents the circular velocity in each 
ring r, while power in other coefficients represent deviations 
from circular motion. On the other hand, the zeroth-order 
term, Ag gives the systemic velocity of each ring. The veloc- 
ity dispersion field is an even moment {^even) of the veloc- 
ity distribution, so that its kinematic analysis is identical 
to traditional surface photometry. In an ideal rotating disk 
the velocity dispersion has to be constant along each ring 



4.2. Kinemetry of LlRGs sample 

In our analysis, the position of the galaxy's center is con- 
sidered at the peak of the Ha flux intensity map. In the two 
class objects, the Ha peak is in good positional agreement 
with the center of symmetry of the velocity field (i.e. kine- 
matic center). For the class 2 objects, the kinematic center 
of the velocity field is not so well defined, but still the Ha 
peak is in a reasonable symmetric position as well. The po- 
sition angle F is left free to vary. The flattening q has also 
been left to vary, but within a physically meaningful range 
(i.e. 0.2-1). This allows us to consider general cases, such as 
tilted/wrapped disks. In Sec. 4.3.1. (i.e., Fig.lTU]) we discuss 
further the effects of the different choices regarding these 
input parameters. 

On figures [5] - [8] we present for the four galaxies the 
original kinematic maps (i.e., velocity field and velocity 
dispersion) along with their reconstructed maps (i.e., ob- 
tained using all the coefficients of the harmonic expansion 
up to the fifth order corrections). The residual maps (i.e., 
original - model) are also shown with a typical RMS of 10 
km/s. Therefore, in all four cases, the fits are good and the 
reconstructed maps recover the properties of the original 
data with great detail. W e also present the behavior of the 
kinematic parameters defined in the previous section (i.e., 
position angle, F, and the fiattening, q) along with some 
of the kinemetry coefficients (ki and ks normalized to the 
ki). 

As for the ellipse parameters, in general we find a good 
optimization for F at the different radii for the four objects. 
As for q, there are several cases which reach the boundary 
of the range of acceptable values. In general this happens 
in the innermost regions where the rotation component, 
and therefore the amplitude of the sinusoidal velocity pro- 
file along the ellipse, are small. Therefore, the irregularities 
in the velocity field are relatively more important when 
compared with the (rotational) amplitude of the velocity 
profile. Since the problem of finding the best ellipse (i.e., 
best q and F combination) can be quite degenerate with 
many local minima, a low q (high ellipticity) facilitates the 
minimization of X^ since the latter is sensitive to a change 
in F. On the contrary, if q is high (circular ellipse), the fit 
(x^) is insensitive to a change in F making more difficult 
the x^ minimization process. As a result, q tends to have 
low values in the inner regions, reaching in some cases the 
minimum acceptable value. In addition the q parameter is 
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less constrained at low radii as a consequence of the rela- 
tively fewer data points involved and the seeing smearing, 
showing larger errortQ. This does not mean, however, that 
the harmonic expansion to describe the kinematics is not 
well constrained. Indeed, the reconstructed velocity maps 
are in excellent agreement with the real velocity field for 
the inner regions (i.e., low q values), even for the inner- 
most region (see figures [5] -[5]). For large r, the amplitude of 
the velocity profile along the ellipse is larger, and more data 
point contribute to constrain the ellipse, being less sensitive 
to local kinematic irregularities and, therefore, providing in 
general a better optimization of the parameters. At large 
radii we also find good agreement between the fit and the 
data. 

It follows some comments on the kinemetry results for 
each galaxy of the sample. 

IRAS F11255-4120 

The position angle F is quite stable (~ 200° ± 15) over 
most the sampled radius. In the inner part (i.e. r < 5"), 
q has values of 0.2-0.4 and it has some instabilities at 
r ~ 7", likely due to the inner bar structure. However 
the expansion recovers very well the data in all this re- 
gion, ki is mainly dominated by rotation (i.e, Bi) and 
increases radially up to a value of 140 km/s. The ks/ki 
term, which measures the small scale kinematic asym- 
metries, has low values up to a maximum of 0.1, point- 
ing some minor peculiarities for r < 7", likely due to the 
bar structure. These asymmetries become smaller in the 
outer part. The error bars are not visible being smaller 
than the black points in the plot. 

IRAS F10567-4310 

The position angle F is quite stable over the sampled 
radii, with values between 210 and 250 degrees. The q 
parameter reaches the value of 0.2 in the inner region 
(i.e., r < 4"), where the reconstructed map shows a very 
good agreement with the data, ki is mainly dominated 
by rotation (i.e, Bi) and increases radially up to a value 
of 140 km/s, revealing a possible minor anomaly at r ~ 
5". The ks/ki term has very low values, with a maxi- 
mum of 0.05 over the sampled radii, meaning that this 
object is close to an ideal rotating disk structure. 





Orig Vel 





fV^- 





^-^fe) 





'^'^HH 



IRAS 11255-4120 
Rec Vel 



Residuals Vel 




-10 10 
arcsec 
Orig Vel Disp 



-10 10 
arcsec 
Rec Vel Disp 



-10 10 
arcsec 

Residuals Vel Disp 




-10 10 -10 10 -10 10 

arcsec arcsec arcsec 

IRAS Fl 1255-41 20 ( Q, F free, C0VER=0.7, 10 poroms) 



300 

' 250 
■ 200 



0.8 
0.6 
0.4 
0.2 
0.0 


r 




^ 


J\j 

• 


•-•--_,-•- --* 


} • . 


H 






















„l20 
..^ 100 
E 80 
^ 60 
- 40 


L 


,-*-;' 


f^~ 








^ 


















0-10 

0.05 

0-00 

-0-05 


r 




,■-- mr 


-«--l 

'v^^- 


-..^^.,-.-. 


* . . 


^ 



R[orc 



Fig. 5. Upper panel: Maps of Ha Gaussian fit velocities (top 
left). Ha Gaussian fit dispersion (bottom left) and their respec- 
tive reconstructed (middle) and residual (data-model) maps (on 
the right) for IRAS F11255-4120. Lower panel: Radial profiles 
of the kinematic properties, obtained using kinemetry program. 
The position angle F and the fiattening q of the best fitting el- 
lipses as well as the first and the fifth order Fourier terms (re- 
spectively, ki and ks) are plotted as a function of the radius. 



IRAS F04315-0840 

The position angle F typically spans from 350° up to 
380if|, except for (r ~ 3") which drops to around 320° 
as a consequence of one of the two distinct redshifted 
peaks present in the velocity field map. For most of the 
radii q is close to 1 (circle-ellipses), ki increases radially 
up to a value of 120 km/s (r < 5") and then decreases to 
reach a value of 100 km/s at r ~ 15". The ks/ki term 
is between and 0.3, with an constant increasing be- 



^ Note that kinemetry provides the errors on the ellipse pa- 
rameters (i.e., q, F), which are determined by the Levenberg- 
Marquardt least-squares minimization (MPFIT) fit with the 
formal la uncertainties computed from the covariance matrix. 
Similarly, the harmonic terms have their formal la errors esti- 
mated from the diagonal elements of the corresponding covari- 
ance matrix obtained with a linear least-square fit. 

^ The position angle is measured from the North (0°=360°) 
anti-clockwise. In order to make the two values more easily com- 
parable we choose to add 360° in one case. 



haviour. This illustrates that departures from rotation 
are mainly found in the outer parts. 

IRAS F21453-3511 

The position angle F almost keeps constant to ~ 240° 
for r > 6" while it shows lower and more unstable values 
for r < 6". The q reaches 0.2 in the region inside a 
radius of r < 6", mainly motivated by the relatively 
small rotation component in this region. However the 
agreement between the kinematic fit and data is very 
good for these radii. At r =7" q reaches a maximum ('^ 
1) to decrease constantly to a value of 0.4. ki does not 
reveal any strong rotational component ranging between 
10 up to 40 km/s. It increases in the inner region for r < 
4" and then remains constant. The ks/ki term ranges 
between and 0.2, showing a quite stable trend. For 
radii larger than 10" it shows an opposite behavior with 
respect to the ki term. In this object the deviations from 
pure rotation are not as high as those of IRAS F04315- 
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Fig. 6. Upper panel: Maps of Ha Gaussian fit velocities (top 
left), Ha Gaussian fit dispersion (bottom left) and their respec- 
tive reconstructed (middle) and residual (data-model) maps (on 
the right) for IRAS F10567-4310. Lower panel: Radial profiles 
of the kinematic properties, obtained using kinemetry program. 
The position angle F and the flattening q of the best fltting el- 
lipses as well as the flrst and the fifth order Fourier terms (re- 
spectively, ki and ks) are plotted as a function of the radius. 



Fig. 7. Upper panel: Maps of Ha Gaussian fit velocities (top 
left). Ha Gaussian fit dispersion (bottom left) and their respec- 
tive reconstructed (middle) and residual (data-model) maps (on 
the right) for IRAS F04315-0840. Lower panel: Radial profiles 
of the kinematic properties, obtained using kinemetry program. 
The position angle F and the flattening q of the best fltting el- 
lipses as well as the flrst and the flfth order Fourier terms (re- 
spectively, ki and ks) are plotted as a function of the radius. 



0840, but they are higher than those for the class 
objects. 

In general we have found that for class sources the 
higher order deviations (from pure rotation) are small (i.e., 
ks/ki < 0.1) while for class 2 objects the deviations are 
higher, mainly in the outer regions (ks/ki < 0.4). For all of 
them the rotation curve (i.e., ki parameter) seem to charac- 
terize these objects as rotating. As discussed before, these 
galaxies do not show extreme asymmetries in their kine- 
matic maps and this is confirmed by the kinemetry results 
obtained so far. 

In the following sections we will study several kinematic 
criteria with the aim of better classifying these systems. 

4.3. Sample of local LIRG systems in the [cja - Va] plane 

In order to reveal the presence of rotational/non-rotational 
motions within the dynamics of the gas in each galaxy we 
will first consider the same criteria as the one proposed by 
SOS. It is worth mentioning that they compare two main 



classes of systems: those which have suffered a recent ma- 
jor merger event (i.e., mergers) and those without sings of 
interacting or merger activity (i.e., disks). For further de- 
tails see SOS. They define the asymmetries in the velocity 
and velocity dispersion fields as: 



^asym 



'^avg.v 



' asym 



^avg,a 

Bi^v 



(4) 



where kavg,v = {k2,v + h.v + ki,v + k^.v)/^ and kavg,a = 
(^i,cr + fe,cr + ^3,cr + fc4,(T + fc5,(T)/5. For an ideal rotating disk, 
we expect the velocity profile to be perfectly antisymmetric 
where the Bi term would dominate the Fourier expansion, 
while the velocity dispersion map is expected to be perfectly 
symmetric and therefore all terms except Aq would vanish. 

In Fig. ini we show the results for our galaxies in the 
[aa — Va] plane. As expected, class objects have lower val- 
ues of the asymmetries than class 2 objects. Indeed, looking 
at their velocity and velocity dispersion maps, the class 
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Fig. 8. Upper panel: Maps of Ha Gaussian fit velocities (top 
left), Ha Gaussian fit dispersion (bottom left) and their respec- 
tive reconstructed (middle) and residual (data-model) maps (on 
the right) for IRAS F21453-3511. Lower panel: Radial profiles 
of the kinematic properties, obtained using kinemetry program. 
The position angle F and the flattening q of the best fltting el- 
lipses as well as the flrst and the fifth order Fourier terms (re- 
spectively, ki and ks) are plotted as a function of the radius. 



maps (Figs. [5] - [S]) resemble to those of an ideal rotating 
disk (i.e., 'spider diagram' structure for the velocity field 
and centrally peaked velocity dispersion map) vifhile class 
2 objects present more distorted velocity fields and irreg- 
ular dispersion maps (Figs. [7]- 15]). Therefore, taking into 
account that the objects in the present sample were classi- 
fied as class or class 2 on the basis of pure morphological 
arguments, we can conclude that their morphological and 
kinematics classification in the [ctq - v^] plane are consis- 
tent. 

In order to analyze the robustness of the results in Fig. 
[51 we will analyze in the following sections their dependence 
on the input parameters considered as well as in the error in 
the radial velocities and velocity dispersion measurements. 
Note that kinemetry uses the full 2D kinematic information 
of the velocity field and velocity dispersion map, allowing 
us a good characterization of the asymmetries, something 
which is more efficient than using ID parameters such as 
Yc/o'c and Y shear /^ as sccn in section 3. 
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Fig. 9. Asymmetry measure of the velocity Vasym and veloc- 
ity dispersion cTasym fields for the four galaxies analyzed here. 
Symbol types distinguish among difi'erent systems: the square 
represents IRAS F10567-4310, the triangle is IRAS F11255- 
4120, the circle is IRAS F04315-0840 and diamond represents 
IRAS F21453-3511. Blue and green symbols stand for difi'erent 
morphological types, respectively class and class 2 objects. 



4.3.1. Dependence on the input parameters 

As pointed out by K06, a good choice of input parameters 
is important in order to avoid an artificial overestimation 
of the asymmetries. To perform a kinemetry analysis we 
require to set the dynamical center, as well as to specify 
different levels of constrains of the input parameters. For 
instance, the kinematic position angle F and the flattening 
q of the ellipses can be fixed or left free to vary for the fitting 
at different radii. The COVER parameter, which controls 
the radius at which the process stops by settingthe fraction 
of the ellipse that has to be covered by dat^j, has to be 
assigned too. 

In Fig. [TU]we present the results in the [aa - Va] plane 
for our galaxies for different sets of input parameters. In 
each panel one parameter at a time is changed: respec- 
tively, the COVER (top- left), the position angle F (or PA, 
top-right), the flattening q (bottom- left) and the CENTER 
(bottom-right) of the ellipses. We can find that the results 
are stable to a reasonable choice of values. Indeed, looking 
at the COVER panel, the four galaxies give similar results 
up to cover=0.5; IRAS F11255-4120 deviates significantly 
for COVER=0.3. IRAS F21453-3511 seems to be sensitive 
to the galaxy center but it maintains within the region of 
high-asymmetries in all the (extreme) cases considered. As 
shown, the choice of free or fixed position angle F or fiat- 
tening q do not affect so much the final results. In par- 



" For instance, COVER = 0.7 means that if less than 70% of 
the points along an ellipse are not covered by data the program 
stops. This value makes sure that kinematic coefficients are ro- 
bust. Sometimes it is necessary to relax this con dition especially 
when reconstructing maps. See further details in lKrainovic et ahl 
(120061). 
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4.3.2. Monte Carlo simulations 
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Fig. 10. Results for the asymmetry measures when different 
sets of input parameters are considered with 10 harmonic terms 
analysis. Symbol types and colors distinguish among the differ- 
ent systems as explain in the legend. In each panel one parameter 
at a time is changed: respectively, COVER, F (or PA), q and the 
CENTER of the ellipses. On the top-left: Different results for the 
COVER parameter (i.e., 0.3, 0.5, 0.7, 0.9) where the size of the 
symbol is proportional to the COVER value (the biggest sym- 
bol corresponds to the highest cover value and vice versa). On 
the top-right: The results obtained considering a position angle 
r free to vary (filled symbols) or fixed to its mean value (empty 
symbols) are shown. On the bottom-left: The results shown are 
obtained when considering q free to vary and constant with the 
radius to its mean value. On the bottom-right: Results achieved 
choosing 5 different centers for each galaxy: filled symbols rep- 
resent results obtained from the 'standard' analysis (Ha flux 
peak), empty symbols the results from shifting their center of 1 
pixel (horizontally and vertically, with respect to the Ha peak 
pixel corresponding to a shift of 0.95"). 



ticular the computed asymmetries are quite insensitive to 
the choice for q, especially for the class galaxies. In gen- 
eral, the results obtained with a fixed T / q are somewhat 
higher due to the fact that there is less degree of freedom 
such that the asymmetries cannot be well accounted by the 
fitting. Obviously the proper choice of free / fixed parame- 
ters depends on the level of S/N (i.e., if the deviations are 
associated to a true feature or noise). 

According to these results we have selected the follow- 
ing set of input parameters for the remainder analysis (i.e., 
COVER — 0.7, r is completely free to vary and q free to 
vary in the range [0.2 — 1], center = Ha flux peak). 



In order to analyze the dependence of the kinemetry re- 
sults to the uncertainties in the radial velocities and velocity 
dispersion values, we measure the probability distribution 
function (PDFs) of the asymmetries in these systems using 
Monte Carlo (MC) simulations, as done in SOS, since the 
kinemetry method does not lend itself to a straightforward 
error propagation. 

Therefore, for each template, we create 150 different re- 
alizations of the moment maps (i.e., velocity field and ve- 
locity dispersion) based on their corresponding error maps. 
These error maps correspond to the measurement errors of 
the velocity moments, as derived when fitting the kinemat- 
ics from the data cube along with wavelength calibration 
errors, as explained in section 2.4. For each moment map 
we perturb the observed data points by randomizing them, 
using Gaussian noise parameterized by the measured Icr er- 
rors. The new maps created are then used to rerun kineme- 
try and apply the analysis described before. In figurefTTjthe 
results are shown. They show that the results in the [a a - 
Va] plane are stable to the velocity errors, with relatively 
well defined regions for disks and merger galaxies. 

4.4. Total kinematic asymmetry of disk, on-going and 
post-coalescence merger systems 

As done in SOS we compute the total kinematic asymme- 
try Ktot border in order to separate class (disks) and 
class 2 (post-coalescence merger) objects. For the four ob- 
jects we obtained a mean value of < Ktot >= 0.135: this 
is a significant lower value than that obtained in SOS (i.e., 
< Ktot >= 0.5). The difference could arise from the fact 
that our mergers are in a post-coalescence phase showing 
relaxation in the innermost regions with a large rotation 
component, while in SOS they are mostly pre-coalescence 
merger systems dominated by dispersion. Also, the limit de- 
fined in SOS comes out after excl uding IRAS 12112+0305, a 
pre-coalescence merger pair (i.e.. lGarcia-Marin et al.ll2009() 
that looks like a disk at high redshift. If the classification 
of this object is considered, the total kinematic asymmetry 
border derived for the whole SOS sample would have been 
considerably lower (i.e., < Ktot >^ 0.3, see Figure 5 in SOS) 
and therefore the discrepancy with our finding reduced. 

In any case, the fact that the border disk/merger for 
the SOS sample is so dependent on the classification of a 
single object (i.e., IRAS 12112+0305) illustrates its rela- 
tively large associated uncertainties. A reduction of the to- 
tal asymmetry border Ktot, as our results (and SOS, after 
reclassifying IRAS 12112-1-0305) suggests, would imply that 
the relative frequency oi disks and mergers in a given sam- 
ple changes, increasing the fraction of mergers. Changing 
the relative frequencies of disk/merger have obvious im- 
plications when interpreting the data in terms of the dif- 
ferent evolutionary scenario mentioned in the introduction. 
On the other hand our sample is admittedly too small and 
formed by objects with relatively homogeneous properties 
(i.e., LIRGs classified as disks and post-coalescence merg- 
ers) for a robust determination of general use. Therefore, 
further efforts to constrain its value as well as to understand 
how it depends on several instrumental and observational 
factors are required. 
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Fig. 11. On the top: Asymmetry measures of the velocity 
Vasym and velocity dispersion aasym fields as derived from the 
Monte Carlo realizations for the four objects. For each source 
150 MC simulations are run but here only the ±3(t results are 
shown. The solid red line indicates the division between class 
(disks) and class 2 (post-coalescence mergers) at < Ktot > = 
0.135. On the bottom: The probability distribution function 
(PDFs) as derived from MC realizations. The empirical delin- 
eation where Kum = 0.135 can cleanly separate the two classes 
(i.e., Kiim area). The symbols are the same as the ones used 
before. 



4.5. A new criterion to distinguisti kinematic asymmetries 
between post-coalescence mergers from disks 

In order to better assess the presence of asymmetries in the 
kinematic maps we will explore a new kinematic criterion. 

As described in iKronberger et all (|2007t ). when consid- 
ering recent (i.e., < 100 Myr after the first encounter) or on- 
going major mergers of equal mass galaxies (i.e.. Milky Way 



type), the inner regionqj of the galaxy are usually affected 
by more chaotic motions, as revealed by a quite irregular 
rotational curve and higher order deviations (i.e., ks/ki) at 
small radii. As the major merger evolves, the inner regions 
rapidly relax into a rotating disk, while the outer parts are 
still out of equilibrium. This implies that the velocity field 
in a post-coalescence system could be dominated by rota- 
tion in the inner regions with large kinematic asymmetries 
in the outer parts, as it is actually observed in our systems. 
Provided that the outer regions retain better the memory 
of a merger event, we propose a criterion which enhances 
the relative importance of the asymmetries at larger radii. 

Indeed, instead of simply averaging the asymmetries 
over all radii (as in SOS), these are weighted according to the 
number of data points used in their determination. Since 
the number of data points of the outer ellipses is larger 
than for the inner ones, the asymmetries found in the outer 
ellipses contribute more significantly to the average for ob- 
taining Yasym and aasym ■ As the number of data points is 
in first approximation proportional to the circumference of 
the ellipse, for practical reasons we used this to weight the 
asymmetries found for the different ellipses. 

The circumferences of the ellipses are computed using 
the truncated 'infinite sum' formula, that is a function of 
the ellipticity (i.e., e = -^/l — q^) and the semi major axis 
of the ellipse (r): 



C(e, r) w 27rr 



1- 



2/ l2-4 



(5) 



The final formula to compute the weighted velocity and 
velocity dispersion asymmetries are respectively: 
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■ ty„ 
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■ L'n 



v^ r 






(6) 
(7) 



where iVis the total number of radii considered, C„ the 
value of the circumference for a given ellipse, the different 
k„ (kj^ and kj^) are the deviations concerning respectively 
the velocity field and velocity dispersion maps, while B^ is 
the rotational curves. We will refer to this approach as the 
'weighted' method and its associated plane W-[cra - v^]). 
In Fig. [12] the results in the W-[CTa - v^] plane for our four 
galaxies are shown where MC have been performed simi- 
larly as above (Sect. 4.3.2). 

The results follow the same general trend but the two 
classes are separated somewhat better than in the un- 
weighted case (Sec. 4.3). Indeed, for class objects the 
weighted velocity asymmetries are lower while, for class 2 
objects are somewhat higher than in the [aa - Vq] plane. 
Therefore, it enhances the fact that post-coalescence merg- 
ers have larger deviations at larger radii with respect to 
pure rotational motions while disks have still lower devia- 
tions than those obtained using [aa - v^]. In this case, the 



^ If a galaxy covers a FoV of 30 kpc x 30 kpc, the regions more 
affected by chaotic motions are those at galactocentric distances 
smaller than 10 kpc. 
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total kinematic asymmetry border, which distinguishes the 
two disks and the two post-coalescence mergers, is charac- 
terized by a mean value of 0.146. 
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Fig. 12. On the top: Weighted asymmetry measures of the 

velocity Va^ym and velocity dispersion a^lm fields as derived 
from the Monte Carlo realizations for the four objects. For each 
source 150 MC simulations are run but only the ±3(7 results are 
shown. The solid red line indicates the empirical division be- 
tween disks and post-coalescence mergers at < Ktot > ~ 0.146. 
On the bottom: The probability distribution function (PDFs) 
as derived from MC realizations. The empirical delineation Kum 
= 0.146 can well separate the two classes along with a large range 
of other values (i.e., dashed Kum area). 



4.6. Angular resolution / Redshift dependence 

The angular resolution effects on the distortions of the ve- 
locity fields produced b y mergers have been discussed by 
iKronberger et ahl (|2007r ) on the basis of simulated velocity 



fields as a function of redshift (i.e., < z < 1). They found 
that for large (Milky Way type) galaxies the distortions are 
still visible at intermediate redshifts but partially smeared 
out, while for small galaxies even strong distortions are not 
visible in the velocity field at z ?» 0.5. 

iGonqalves et al.l ()2010t ). simulating LB As at redshift 
z ^ 2, found that, in general, galaxies at high redshift 
present smaller values of Kasym, i-e., they appear more 
disky than they actually are. The percentage of galaxies 
classified as mergers drop from ~ 70 % to ^ 38 % from 
z = to z = 3 according to their simulations. 

In order to investigate the resolution effects on our re- 
sults we simulate to 'observe' these systems at z = 3 with 
a typical pixel scale of 0.1" (the same pixel scale of the 
IFU NIRSpec/JWST). At this redshift the current FoV of 
our images is about 1" x 1" with a typical scale of 7.83 
kpc/arcsec, assuming a ADCM cosmology with Hq — 70 
km/s/Mpc, r^M = 0.3 and JIa = 0.7. The simulated maps 
are shown in Fig. 1131 We apply kinemetry using these maps 
and, following the same procedures as before, obtain the 
results shown in the [aa - Va] (Fig. [H|) and W-[(Ta - Vq] 
(Fig. [15]) planes. The [aa - v^] plane shows the expected 
trend, where both classes appear more symmetric when ob- 
served at high redshift. In this case a lower value of the total 
kinematic asymmetry }iasym border is derived (red dashed 
line), as expected (i.e., Kasy7n = 0.096). Thus, shifting the 
sample from z = to z = 3, the frontier between disks / 
post-coalescence mergers in the [aa - Va] plane changes from 
0.135 to 0.096. 

The W-[aa - Va] plane is less sensitive to resolution ef- 
fects after redshifting our sample at z = 3 (see Fig. [151). This 
is due to the fact that the larger / outer regions are the ones 
less affected by resolution effects. As with this criterion the 
associated asymmetries weight more than those present at 
inner radii, the computed asymmetries (i.e., Vasym, cTasym) 
are less affected by resolution. Therefore, the total kine- 
matic asymmetry distinguishing the two classes changes 
only from 0.146 to 0.130 between z = and z = 3. 

Summarizing, resolution effects tend to smooth kine- 
matic deviations making objects to appear more disky than 
they actually are. This effects are more significant when 
analyzing the kinematic asymmetries in the (unweighted) 
[o'a -Va] plane than in the W-[(7q - Vq] one. In particular, 
when compared our local seeing limited observations with 
simulated data at z = 3 (and 0.1"/spaxel), the total kine- 
matic asymmetry border value is reduced by a 30% from 
z = to z = 3 in the [aa - Va] plane, while it is only shifted 
by 11% in the W-[aa - Va] plane. 

Comparing our results with those obtained in S08 when 
considering local spirals and toy-disk models as observed 
at high redshift, we notice that the low value of our K^tot 
(~ 0.1) classifies those galaxies as mergers and, on the 
other hand, the frontier defined in S08 classifies our post- 
coalescence mergers and IRAS F12112-I-0305 as disks. This 
illustrates that the definition of merger is a crucial point 
to define the asymmetry frontier. Indeed, a lower value of 
the fraction merger/disk (i.e., a higher value of the Ktot 
parameter) can be derived if only pre-coalescence on-going 
mergers are considered as 'true' mergers. The weighted cri- 
terion proposed here should be applied to larger and more 
diverse samples in order to understand the uncertainties 
associated to this type of classifications. 
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Fig. 13. Ha maps for the four objects as observed at z = 3 
with a spatial samphng of 0.1". The flux intensity, velocity 
dispersion a (km/s) and velocity fields v maps (km/s) for 
the main component are shown. The flux intensity maps 
are represented in logarithmic scale and in arbitrary flux 
units. All the images are centered using the Ha peak and 
the iso-contours of the Ha flux are over-plotted. 



5. Conclusions 

This paper presents the results from spatially resolved kine- 
matics of four local (z ~ 0.016) luminous infrared galax- 
ies (i.e., LIRGs) observed with the VLT/VIMOS IFU as 
part of a larger project to characterize the properties of 
(U)LIRGs on the basis of optical and infrared Integral Field 
Spectroscopy. The four galaxies are at a similar distance 
(~ 70 Mpc) and, on the basis of its morphology, two of 
them have been classified as isolated disks and two as post- 
coalescence merger objects (see Paper HI). The velocity 
field and velocity dispersion maps are analyzed with the 
aim of studying in detail their kinematics. The kinemetry 
method (developed by Krajnovic and coworkers) is used 
to characterize the kinematic asymmetries of these galax- 
ies and several criteria are discussed to distinguishing their 
status. We can draw the following conclusions from this 
study: 

— The general kinematic properties of the four LIRGs are 
consistent with their morphological classification: iso- 
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Fig. 14. The comparison between local and 2 = 3 results in 
the [aa - Va] plane. Empty symbols represent low-z results while 
filled ones are for 2 = 3. Symbols have the same meaning as in 
Fig. [5] The red dashed line is the high-z frontier with a value 
of 0.096 as explained in the text while the dotted one is for the 
local analysis (i.e., 0.135). 
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Fig. 15. Comparison between local and 2 = 3 results in the 
W-[(Ta - Va] plane. Empty symbols represent low-z results while 
filled ones are for z=3. The red dashed line represents the high-z 
frontier, while the dotted one is for the local case. 



lated disks reveal quite regular velocity fields and cen- 
trally peaked velocity dispersion maps consistent with 
a single rotating disk interpretation while the remain- 
ing two show departures from this behavior. In partic- 
ular we found that post-coalescence mergers show more 
irregular velocity fields and velocity dispersion maps 
showing off-nuclear dispersion peaks (up to 220 km/s 
at 2.4 kpc from the H„ peak for IRAS F04315-0840) or 
nuclear asymmetric structures; 

We have found double peaked emission line profiles in 
the inner regions of the four galaxies. The secondary 



15 



Bellocchi et al.: Kinematic asymmetries of disks and post-coalescence mergers 



broad components (i.e., ct ~ 70 - 450 kni/s) are in 
all cases blue-shifted (Aw ~ 50 - 150 kni/s) and, tak- 
ing into account the large velocities involved, they can 
be well explained by the presence of an outflow in a 
dusty environment. The pattern of the 2D kinematic 
maps of the secondary broad component for the two 
post-coalescence mergers (with kinematic axes perpen- 
dicular to those of the main component) further sup- 
ports this interpretation. In the particular case of IRAS 
F043 15-0840 the broad component is found over a quite 
extended area (^ 2.4 kpc x 2.7 kpc); 
The Vc/o'c parameter classifies our sources as rotation 
dominated. Similar results are obtained when using the 
quantity y shear l'^- This shows that our post-coalescence 
mergers have a large rotation component and the sole 
use of these parameters does not allow us to discriminate 
their kinematic differences with respect to disks; 
When the full 2D kinematics information provided by 
the spatially resolved velocity field and velocity disper- 
sion maps is considered, the kinematic asymmetries are 
well characterized with kinemetry, making the morpho- 
logical and kinematic classification consistent for the 
four objects. Disks have lower kinematic asymmetries 
than those derived for post-coalescence mergers; 
We have explored a new criterion to characterize the 
kinematic asymmetries using kinemetry. In particu- 
lar, we introduce a new weighting method which gives 
weight to the kinematics of the outer regions when com- 
puting the total asymmetries Vasym and aasym- This 
is motivated by the fact that post-coalescence mergers 
show relatively small kinematic asymmetries in the in- 
ner parts as a consequence of the rapid relaxation into 
a rotating disk, with the outer parts being still out of 
equilibrium (i.e., larger asymmetries). The fact that the 
'frontier' between disks and post-coalescence systems 
only changes by 11% when considering the local and 
high-z cases suggests that this new criterion is less de- 
pendent on angular resolution effects. Thus, the W-[(7a 
- Va] plane differentiate in a more robust way disks and 
post-coalescence mergers. 

Classifying disk/merger systems using kinemetry is 
difficult, and it obviously depends on the definition of 
merger. The 'asymmetry frontier' strongly depends on 
the 'type' of mergers considered: if only pre-coalescence 
on-going mergers are considered as 'true' mergers a 
lower value of the fraction merger/disk systems can 
be derived. The weighted criterion proposed here helps 
to characterize in a more robust way post-coalescence 
asymmetries. 

Using previous criteria to classify disks/mergers, our 
two post-coalescence systems would have been classi- 
fied as disks. This suggests that the ratio of mergers to 
disks systems at high-z may have been underestimated. 
Larger and more diverse samples are required to confirm 
this conclusion. 
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Table 3. Kinematic parameters for the sample. 



Galaxy ID 
IRAS code 


Type of fitting ^ 


■ b 

degree 


(km/s) 


(fem/s) 


Vc/o-c 


{kraj s) 


(km/s) 


Veh = a./S 


(kpc) 


M^„„ 
(Mo) 




F11255-4120 
F11255-4120 


Ic 

2c 


53° ± 20° 
53° ± 20° 


187 ± 72 

188 ± 76 


104 ±3 
83 ±2 


1.8 ± 0.7 
2.3 ± 1.0 


125 


47 


2.7 


2.40 ±0.57 


(2.9 ±0.8) • 


ioi« 


F10567-4310 
F10567-4310 


Ic 
2c 


36° ± 4° 
36° ± 4° 


253 ± 32 
255 ± 30 


56 ±2 
55.0 ± 1.4 


4.5 ± 0.8 

4.6 ± 0.7 


120 


41 


3 


3.20 ±0.40 


(1.7 ±0.3) • 


10i« 


F04315-0840 
F04315-0840 


Ic 

2c 


n.c. 
n.c. 


157 ±3 
162 ±5 


110 ±5 
69± 3 


1.4 ± 0.1 
2.3 ± 0.2 


51 


51 


~ 1 


0.94 ±0.14 


(7.9 ± 1.9) 


■10» 


F04315-0840 
F04315-0840 


Ic 
2c 


29° ± 3° 
29° ± 3° 


324 ± 40 
335 ± 45 


110 ±5 
69 ±3 


2.9 ± 0.5 
4.8 ± 0.9 


51 


51 


~ 1 


0.94 ±0.14 


(7.9 ± 1.9) 


•lO^ 


F21453-3511 
F21453-3511 


Ic 
2c 


n.c. 
n.c. 


90 ± 1 

78 ±2 


69 ±3 
61 ±1 


1.30 ± 0.07 
1.28 ± 0.05 


132 


65 


2 


2.97 ±0.82 


(1.9 ±0.6) ■ 


10" 


F21453-3511 
F21453-3511 


Ic 
2c 


24° ± 2° 
24° ± 2° 


222 ±21 
192 ±21 


69 ±3 
61 ± 1 


3.2 ± 0.5 
3.2 ± 0.4 


132 


65 


2 


2.97 ±0.82 


(1.9 ±0.6) ■ 


10" 



d 



W 



O 



B 



a- 
o 



Notes. " 1 component fitting (Ic) and 2 component fitting (2c, referees to tlie main (systemic) component of tfie two-Gaussian fit) is consider 

for tfie four galaxies. Inclination of the galaxy; n.c. is when no inclination correction is applied. '^ Circular velocity derived as the half of the o 



o 



observed peak-to-peak velocity from the Ha kinematic. Corrected and no corrected values for the inclination are shown. Central velocity dis- £. 

persion as derived from the ana maps. "^ Velocity shear not corrected for the inclination of the galaxy. •' Global velocity dispersion in the whole galaxy. 
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Appendix A: Other kinemetry-based criteria for 
detecting asymmetries 

In this appendix we describe other criteria that could be 
useful to distinguish different systems. Indeed, other au- 
thors tried to find out criteria able to kinematicahy classify 
galaxies at interme diate redshift (i.e., z ~ 0.6): for example, 
iFlores et al.l ()2006l ) developed a simple kinematic classifica- 
tion scheme for distant galaxies based on their 3D kinemat- 
ics and their morphology in the ACS F775W images. They 
model the velocity field of a rotating disk that matches the 
observed velocity gradient to obtain the expected cr-map 
corresponding to the observed velocity field and then com- 
pare the observed and model cr-maps, estimating whether 
the observed kinematics are consistent or not with a rotat- 
ing disk (i.e., a rotating disk should show a well defined 
peak in the center of the a map). The spatial separation 
(Ar, in pixel) between the peaks in the two cr-maps and 
the relative difference (e) between the amplitudes of the 
modeled and observed a peaks are considered to classify 
the objects. They found that rotating disks have locations 
near Ar ~ e '-~^ 0, while galaxies with anomalous velocity 
fie lds are clear l y off set. The same criterion has been applied 
in I Yang et ^11200 8). 

Here we propose some criteria, based again on the 
kinemetry method: 

1. In this case we consider the normalized higher order de- 
viations <k5 / ki>„ of the velocity field with respect to 
those of the velocity dispersion map Oasym , since the ks 
term indicate s complex structure in the velocity map, 
as noticed in iKrainovic et ahl (J2006l) . Small ks/ki val- 
ues should be found in disk-like structures while high 
ks/ki are expected for mergers, where more complex 
structures are supposed to exists. This is confirmed by 
our results: class objects have <k5 / ki> close to zero, 
the amplitude of fci is substantial, the position angle T 
and fiattening q remain quite constant such that they 
can be classified as disk-like objects. On the other hand, 
class 2 objects show a rise of the <k5 / ki> term up 
to 0.6 at outer radii, where peculiar structures can be 
identified. In Fig. lA. li the results obtained from apply- 
ing MC simulations (150 for each object). This criterion 
allows to distinguish the class / 2 in similar way as in 
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the [(Ta - Va] plane, but somewhat worst than in W-[cr£i 
- Va] plane. 

3ct results 



0.10 



A 

> 

V 



0.01 




0.01 



3ct results 



0.10 



(kSI) 



Disks 

P.c. Mergers 

Ki:m area 



Q 
CL 




0.0 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 

Fig. A.l. On the top: Asymmetry measures of the <k5/ki > 
term and velocity dispersion al^lrn. as derived from the Monte 
Carlo realizations for the four objects. For each source 150 MC 
simulations are run but the ±3(t results are shown. The solid red 
line indicates the empirical division between disks and mergers 
at <Kiim > = 0.085. On the bottom: The probability distri- 
bution function (PDFs) as derived from the respective MC real- 
izations. The empirical delineation well separate the two classes 
defining a small range of possible values (i.e., Y^um area). Small 
deviations <k5/ki > are found for disks while higher values for 
mergers, consistent with our expectation. 



map and the mean values of the -Bi,i, term over all the 
radii. The next formula describes the two asymmetries: 



^asym 



RMS 



■' asym 



Bi^v 



where the RMS gives an estimate of the deviations of 
the system from the ideal case of simple rotation in the 
systerrO- In Fig- IA.2I the results in 5 - [aa - Va] plane 
are shown. Comparing this plot with that shown in Fig. 
[21 we notice that the RMS/<Bi> term tends to be 
higher than the Vasym values while the o'^asym'' range 
smaller than the (y^asynT''^ range. The o'^asym'^ does not 
seem to characterize well the kinematic asymmetries of 
these systems (i.e., quite distorted for the IRAS F04315- 
0840). On the other hand, the RMS/<Bi> term by it- 
self can distinguish the two classes relatively well, where 
a border dividing them can be estimated with a value 
of ~ 0.2. 



"bo results 



(5 terms) 




0.10 



'osym 



Fig. A. 2. Subtracting the circular motion (i.e., Bi coefficient) 
from the data we end up with the velocity residual maps, char- 
acterized by mean residuals close to zero and RMS lower than 
20 km/s for the whole sample. Having the four objects almost 
the same RMS, this plot clearly implies higher circular velocities 
for the class galaxies while lower Bi term is derived for the 
class 2 ones. Their rotation is well shown in Figs. ((5]-|8l). 



In some cases the data quality does not allow us to 
apply sophisticated criteria to look for asymmetries. In 
these cases a simple criterion may be the only possible 
option. Therefore, we have studied the potential of the 
Simple-[aa - Vq] plane (hereafter, S-[aa - v^,]), where 
only the first order correction in the harmonic expansion 
is considered and asymmetries are defined as the ratio 
between the root mean square (RMS) of the residual 



Residuals are from measured velocities and Bi coefficient. 
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